Strain effect on the optical conductivity of graphene 
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Within the tight binding approximation, we study the dependence of the electronic band structure 
and of the optical conductivity of a graphene single layer on the modulus and direction of applied 
uniaxial strain. While the Dirac cone approximation, albeit with a deformed cone, is robust for 
sufficiently small strain, band dispersion linearity breaks down along a given direction, corresponding 
to the development of anisotropic massive low-energy excitations. We recover a linear behavior of 
the low-energy density of states, as long as the cone approximation holds, while a band gap opens 
for sufficiently intense strain, for almost all, generic strain directions. This may be interpreted in 
terms of an electronic topological transition, corresponding to a change of topology of the Fermi 
fine, and to the merging of two inequivalent Dirac points as a function of strain. We propose that 
these features may be observed in the frequency dependence of the longitudinal optical conductivity 
in the visible range, as a function of strain modulus and direction, as well as of field orientation. 

PACS numbers: 78.40.Ri, 62.20.-x, 81.05.Uw 



I. INTRODUCTION 

Graphene is an atomic thick single layer of carbon 
atoms in the sp^ hybridization state, which in normal 
conditions crystallizes according to a honeycomb lattice. 
The quite recent realization of sufficiently large graphene 
flakes in the laboratoryi*^ has stimulated an enormous 
outburst of both experimental and theoretical investi- 
gation, due to its remarkable mechanical and electronic 
properties, that make graphene an ideal candidate for 
apphcations in nanoelectronics (see Ref. for a recent, 
comprehensive review). 

Because the honeycomb lattice is composed of two in- 
terpenetrating triangular sublattices, graphene is char- 
acterized by two bands, linearly dispersing at the so- 
called Dirac points. Indeed, this is suggestive of the pos- 
sibility of observing relativistic effects typical of quan- 
tum electrodynamics in such a unique condensed mat- 
ter systent^!^. The presence of massless low-energy ex- 
citations also endows the density of states with a linear 
dependence on energy at the Fermi level, which makes 
graphene a zero-gap semiconductor. This in turn de- 
termines most of the peculiar transport properties of 
graphene, including a minimal, finite conductivity in the 
clean limit at zero temperature^, and a nearly constant 
conductivity over a large frequency interval^'^-. 

Graphene is also notable for its remarkable mechani- 
cal properties. In general, nanostructures based on sp^ 
carbon, such as also nanotubes and fuUerenes, are char- 
acterized by exceptional tensile strengths, despite their 
reduced dimensionality. In particular, recent ah initio 
calculations^ as well as experiments^ have demonstrated 
that graphene single layers can reversibly sustain elas- 
tic deformations as large as 20%. In this context, it has 
been shown that Raman spectroscopy can be used as a 



sensitive tool to determine the strain as well as some 
strain-induced modifications of the electronic and trans- 
port properties of graphen o^°i^^ . 

This opened the question whether applied strain could 
induce substantial modifications of the band structure 
of graphene, such as the opening of a gap at the Fermi 
level, thereby triggering a quantum phase transition from 
a semimetal to a semiconductor. While earlier ah initio 
calculations were suggestive of a gap opening for arbi- 
trary strain modulus and direction^^, both tight-binding 
modelsi^ as well as more accurate ah initio calculationsi^ 
point towards the conclusion that the strain-induced 
opening of a band gap in fact depends critically on the 
direction of strain. On one hand, no gap opens when 
strain is applied in the armchair direction, whereas on 
the other hand a sizeable strain modulus is required in 
order to obtain a nonzero band gap, for strain applied 
along a generic direction. This result is in some sense 
consistent with the overall conviction that the electron 
quantum liquid state in graphene, characterized by low- 
energy massless excitations, lies indeed within some sort 
of 'quantum protectorate', i.e. it is stable against suffi- 
ciently small, non-accidental perturbations^^. 

In this paper, we will be concerned on the effects in- 
duced by applied strain on the optical conductivity of 
graphene. Among the various peculiar transport proper- 
ties of graphene, the conductivity as a function of fre- 
quency ljj and wavevector k in graphene has received 
considerable attention in the past (see Ref. for a re- 
view). The optical conductivity has been derived within 
the Dirac-cone approximation'', and within a more accu- 
rate tight-binding approximation also for frequencies in 
the visible range®. The effect of disorder has been consid- 
ered by Peres et ali^, and that of finite temperature by 
Falkovsky and Varlamov— >. These studies are consistent 
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with the experimentally observed of a nearly constant 
conductivity of ire^ /2h over a relatively broad frequency 
rang e — - ' - P . Such a result demonstrates that impurities 
and phonon effects can be neglected in the visible range 
of frequencies^^. 

Although uniaxial strain will be included in a stan- 
dard, non-interacting model Hamiltonian at the tight- 
binding level, i.e. through the introduction of strain- 
dependent hopping parameterai^, this will nonetheless 
capture the essential consequences of applied strain on 
the band structure of graphene. In particular, the strain- 
induced modification of the band structure at a fixed 
chemical potential may result in an electronic topolog- 
ical transition (ETT)^^ (see Refs. 22,23 for comprehen- 
sive reviews). In metallic systems, this corresponds to 
a change of topology of the Fermi line with respect to 
an external parameter, such as the concentration of im- 
purities or pressure, and is signalled by the appearance 
of singularites in the density of states and other derived 
thermodynamic and transport properties. The effect of 
the proximity to an ETT is usually enhanced in systems 
with reduced dimensionality, as is the case of graphene. 
Here, one of the main consequences of applied of strain 
is that of moving the Dirac points, i.e. the points where 
the band dispersion relations vanish linearly, away from 
the points of highest symmetry in the first Brillouin zone, 
and of deforming their low-energy conical approximation. 
Moreover, for a sufficiently large strain modulus and a for 
a generic strain direction, two inequivalent Dirac points 
may merge, thus resulting in the opening of a band gap 
(at strains larger than the critical one) along one specific 
direction across the degenerate Dirac point. This results 
in a sublinear density of states exactly at the transition, 
which in turn gives rise to an unusual magnetic field de- 
pendence of the Landau levels^!. This may also be de- 
scribed as a quantum phase transition, from a semimetal 
to a semiconductor state, of purely topological origin^, 
characterized by low-energy massless quasiparticles de- 
veloping a finite mass only along a given direction. 

It may be of interest to note that similar effects have 
been predicted also for other low-dimensional systems, 
and that their overall features are generic with respect 
to their detailed crystal structure. In particular, a sim- 
ilar discussion applies to some quasi-two-dimensional 
Bechgaard salts^^i^, as well as to cold atoms in two- 
dimensional optical latticesSAiSS,, which have been pro- 
posed to simulate the behavior of Dirac fermionsSi. 

The paper is organized as follows. In Sec.|TT]we review 
the tight binding model for strained graphene, discuss 
the location of the Dirac points, and the occurrence of 
the ETTs, as a function of strain. In Sec. IIIII we derive 
the Dirac cone approximation for the band dispersions 
in the presence of strain, and discuss the low-energy en- 
ergy dependence of the density of states. This is then 
generalized over the whole bandwidth beyond the cone 
approximation. The formation of band gaps is discussed 
with respect to strain modulus and direction. In Sec. IIVI 
we present the main results of this paper, concerning the 



optical conductivity of strained graphene, and relate the 
occurrence of several singularities in the frequency de- 
pendence thereof to the various ETTs. We summarize 
our conclusions and give directions for future studies in 
Sec. El 



II. MODEL 

Within the tight-binding approximation, the Hamilto- 
nian for the graphene honeycomb lattice can be written 
as 

= ^tfa^(R)6(R-|-5^)+H.c., (1) 

where (R) is a creation operator on the position R of 
the A sublattice, 6(R -I- 5i) is a destruction operator on 
a nearest neighbor (NN) site R -|- 5^ , belonging to the B 
sublattice, and 5^ are the vectors connecting a given site 
to its nearest neighbors, their relaxed (unstrained) com- 
ponents being df^ = a(l,V3)/2, = a(I,-V3)/2, 
Sf^ = a(-l,0), with a = 1.42 A, the equilibrium C-C 
distance in a graphene sheet^. In Eq. ([T]), ti = t{5i)^ 
I = 1,2,3, is the hopping parameter between two NN 
sites. In the absence of strain they reduce to a single 
constant, tg = to, with to — —2.8 eV (Ref. [30l ). 
In terms of the strain tensor^^ 

_ / cos^ 6* - sin^ 6* (1 + i/) cos sin , , 

1^(1 + 1^) cos 61 sine sin^ 6* - j/cos^ 61 y ' ^' 

the deformed lattice distances are related to the relaxed 
ones by 

Se^iI + e)-Sf\ (3) 

In Eq. 9 denotes the angle along which the strain is 
applied, with respect to the x axis in the lattice coor- 
dinate system, e is the strain modulus, and = 0.14 is 
Poisson's ratio, as determined from ab initio calculations 
for graphene^, to be compared with the known experi- 
mental value I' = 0.165 for graphite^. The special values 
6 = and 9 = 7r/6 refer to strain along the armchair and 
zig zag directions, respectively. Fig. [1] shows a schematic 
representation of the strained graphene sheet, along the 
generic direction 9 = 7r/4, for definiteness. 

Eqs. and ^ rely on the assumption that the lat- 
tice structure of graphene responds elastically to applied 
strain, so that the effect of strain on the electronic prop- 
erties can be studied straightforwardly. The robustness 
of the elastic picture is confirmed by recent atomistic 
simulations which have been compared with available 
experiments^. It should be however mentioned that fur- 
ther detailed calculations^i^ii^ show that the formation 
of topological defects under strain, such as Stone- Wales 
defects, can even induce structural phase transitions both 
in graphene and nanotubes. 
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FIG. 1: (Color online) Schematic representation of the 
strained (red) vs unstrained (blue) honeycomb lattice, for 
9 — n/4. Open (closed) symbols refer to the A (B) sub- 
lattices, respectively. 



Let SLe {af ^ ) the strained (unstrained) basis vectors of 

the direct lattice, and he (b^°^ ) the strained (unstrained) 
basis vectors of the reciprocal lattice, respectively, witt^ 



.(0) 



(27r/3a)(l,±\/3). One has = T • bf \ Then, One finds 



from Eq. ([3]), it is straightforward to show that T 
(I + e)"^. Since the wavevectors with and without ap- 
plied strain are connected by such a bijective tranforma- 
tion, we can safely work in the unstrained Brillouin zone 
(IBZ), i.e. k = k(o). 

It is useful to introduce the (complex) structure fac- 
tor in momentum space, as well as the NN hopping and 
overlap functions, respectively defined as 



7k 



e=i 

3 



i=i 

3 



5k 



, See 



(4a) 
(4b) 
(4c) 



Here, the strain-dependent overlap parameters se are a 
generalization of the band asymmetry parameter sq of 
Ref. aO, and are defined as 



(5) 



dr4>{r)4>{r + de) = exp -—^ 



Here, 



(l/V^ag)exp(-rV2a2) is 



normal- 



ized gaussian pseudoatomic wavefunction, with 1/<Jg 



Zg/2^/3a (Refs. IsalSa) , and the value Zg = 11.2 is fixed 
by the condition that the relaxed overlap parameter be 
So = 0.07 (Refs. '_3C|[36|). Correspondingly, the hopping 
parameters tg are defined as the transition amplitudes 
of the single-particle Hamiltonian, Hi = — (fi,^/2m)V'^ — 
Ze^ /r, between two lattice sites being 5i apart from each 
other. Here, Z is chosen so that te = io in the unstrained 



limit. One finds 

2'mag 



Ze^ 



-lo 



se, (6) 



where Io{x) is a modified Bessel function of the first 
kind^^. One finds dte/d6e = 7.6 eV/A for e = 0, which is 
comparable with the value 6.4 eV/A obtained in Ref. [l3 
within Harrison's approacb^^. In the unstrained limit 
(e = 0), Eqs. (gl) reduce to /k tojk and so7k, 
respectively. 

Within the tight binding approximation, the energy 
dispersion relations can be obtained as the solutions i?kA 
of the generalized eigenvalue problem 



-ffkUkA — -EkA'S'kUkA, 



where 



Hi, = 



/k 
/k 

1 .9k 

.9k 1 



-EkA — 



^kT VJ^i^+4Gk|/k| 
2Gk 



(7) 

(8a) 
(8b) 

(9) 



where A = 1 (minus sign) refers to the valence band, and 
A = 2 (plus sign) refers to the conduction band, and 



^ic = 5k/k + .9k/k, 



Gk = l-|5kr- 
The eigenvectors UkA in Eq- can be presented as 



(10a) 
(10b) 



UkA 



cos (pi^x 
-'"'^ sin0kA 



where e'^'' = ./k/|/k|, and 



COS0kA = - 



sm (^kA 



^kAvGk 



l/kP + Gi^El-^ 
l./k| 

|/kP + Gk£;^^^ 



(11) 

(12a) 
(12b) 



with cos((/)k,i — (t>k..2) = 0. In the limit of no strain, one 
finds (/)k,i —> 37r/4 and 0k, 2 ~> 7r/4. Here and below, 
A = 2 when A = 1, and vice versa. 

As already observed in Ref. 36, a nonzero value of the 
overlap parameters se endows the conduction and va- 
lence bands with a finite degree of asymmetry, which 
is here increasing with increasing modulus of applied 
strain and, in general, anisotropic, depending on the di- 
rection of the applied strain. In the unstrained limit 
(e = 0), one recovers the band dispersions of Ref. [s^, 
£;kA ^ ±to|7k|/(l ± so|7k|) (with to < 0). 
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FIG. 2: (Color online) Contour plots of the dispersion rela- 
tions within the IBZ for the valence band, E^i (left panel), 
and conduction band, _Ek2 (right band), Eq. ([9]). Here, we 
are depicting the situation corresponding to a strain modulus 
of e = 0.18 along the generic direction 6 — 7r/4. Solid blue 
lines are separatrix lines and occur at an electronic topologi- 
cal transition, dividing groups of contours belonging to differ- 
ent topologies. Either line passes through one of the critical 
points Ale (^ = 1,2,3), defined as the middle points of the 
IBZ edge (solid black hexagon). 



The band dispersion relations E^x, Eq. are char- 
acterized by Dirac points, i.e. points in k-space around 
which the dispersion is linear, when /k = 0. As a function 
of strain, such a condition is satisfied by two inequivalent 
points ±k£) only when the 'triangular inequalities' 



\tt,-ti,\ < \ti,\ < 



(13) 



are fulfilled^*', with {£i,£2,^3) a permutation of (1,2,3). 
Around such points, the dispersion relations i?kA can be 
approximated by cones, whose constant energy sections 
are ellipses. 

The location of ik^i in the reciprocal lattice satisfies 



cos(k£, • {6i^ -Si^)) = 



- 



(14) 



with (^1,^2,^3) a permutation of (1,2,3). While in 
the unstrained limit the Dirac points are located at 
the vertices of the IBZ (Ref. [j), i.e. k^ — > K = 
(27r/3a, 27r/3V3a) and — k/j — > K' = — K, when either 
of the limiting conditions in Eqs. (jl3p is fulfilled as a 
function of strain, say when tg^ — te-^ + ti^, the would-be 
Dirac points coincide with the middle points of the sides 
of the IBZ, say kz, ^ Mt,. Here, Mi = ^(-l,-\/3), 
M2 = ^(-l,^/3), M3 = §J(1,0). In this limit, the dis- 
persion relations cease to be linear in a specific direction, 
and the cone approximation fails. 

Fig. [2] shows contour plots of E'kA, Eq. at constant 
energy levels. For fixed strain, each of these lines can 
be interpreted as the Fermi line corresponding to a given 
chemical potential. One may observe that the various 
possible Fermi lines can be grouped into four families, ac- 
cording to their topology. In particular, from Fig. [2] one 
may distinguish among (1) closed Fermi lines around ei- 
ther Dirac point ik/^ (and equivalent points in the IBZ), 
(2) closed Fermi lines around both Dirac points, (3) open 



Fermi lines, (4) closed Fermi lines around F = (0, 0). The 
transition between two different topologies takes place 
when the Fermi line touches the boundary of the IBZ 
(solid black hexagon in Fig. and is marked by a sep- 
aratrix line. It can be proved explicitly that the Fermi 
line at the transition touches the boundary of the IBZ 
precisely at either of the hexagon sides midpoints Ali 
{i = 1,2,3), defined above. This situation holds exactly 
also in the presence of overlap {sg 7^ 0), as can be proved 
within group theory^. 

Each separatrix line corresponds to an electronic topo- 
logical transition (ETT)— i^i^^, i.e. a transition between 
two different topologies of the Fermi line. An ETT 
can be induced by several external parameters, such as 
chemical doping, or external pressure, or strain, as in 
the present case. The hallmark of an ETT is provided 
by a kink in the energy dependence of the density of 
states (DOS) of three-dimensional (3D) systems, or by a 
logarithmic cusp (Van Hove singularity) in the DOS of 
two-dimensional (2D) systems. Besides being thoroughly 
studied in metala2^i2^, the proximity to an ETT in quasi- 
2D cuprate superconductors has been recently proposed 
to justify the nonmonotonic dependence of the supercon- 
ducting critical temperature on hole doping and other 
material-dependent parameters, such as the next-nearest 
neighbor to nearest neighbor hopping ratio^, as well as 
several normal state properties, such as the fluctuation- 
induced excess Hall conductivitj*^. In particular, the 
role of epitaxial strain in inducing an ETT in the cuprates 
has been emphasized^. Due to the overall C2V symme- 
try of the underlying 2D lattice in the Cu02 layers, at 
most two (usually degenerate) ETTs can be observed in 
the cuprates. Here, in the case of strained graphene, 
characterized instead by D^h symmetry, we surmise the 
existence of at most three, possibly degenerate, ETTs, 
whose effect on observable quantities may be evidenced 
by the application of sufficiently intense strain along spe- 
cific directions. 



III. DENSITY OF STATES 

Under the conditions given by Eqs. ([13]), the band dis- 
persions, Eqs. ([5]), can be expanded as E^x = i?kA around 
either Dirac point, say k = k/j + q, as: 



P -q-dT V(q-d)2+4Gk„|q-V/J 



where 



2GkD 



(16) 



Eq. (fTS]) defines a cone, whose section E^x = E aX & 
constant energy level E is an ellipse. Its equation can be 
cast in canonical form as 



{qx - qxnY , {qy - qynf 



52 



E' 



(17) 
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FIG. 3: (Color online) Showing the DOS prefactor pi, Eq.[T9] 
normalized with respect to its value pj"' in the absence of 
strain, as a function of the strain modulus e, for various strain 
angles. The strain direction 6 increases from 9 = Q (armchair 
direction, corresponding to the lowest curve) to 6 = n/Q (top- 
most curve). All other cases can be reduced to one of these 
exploiting the symmetry properties of the lattice. 



where the various parameters entering Eq. (|17p are de- 
fined in App. [X] Making use of Eq. p7|) . one can derive 
the low-energy expansion of the density of states (DOS), 
which turns out to be linear in energy. 



p{E)^p,\El 



(18) 



with 



Pi 



-^{tl+tl+tlf-2{t\ + t\ + tt)]-''\ (19) 



where the factor of four takes into account for the spin 
and valley degeneracies. 

Fig. [3] shows the prefactor pi, Eq. as a func- 

tion of the strain modulus e, for various strain angles 
9. One finds in general that pi increases monotonically 
with increasing strain. Such a behavior suggests that ap- 
plied strain may be used to amplify the DOS close to the 
Fermi level. This, in particular, may serve as a route to 
improve known methods to increase the carrier concen- 
tration in doped graphene samples. When the equality 
sign in Eqs. is reached, the prefactor pi in Eq. ([TO]) 
diverges, meaning that the cone approximation breaks 
down. In this case, the band dispersions still vanish, but 
now quadratically along a specific direction through the 
would-be Dirac point, and a nonzero gap in the DOS 
opens around E = Q. 

This behavior is confirmed by the energy dependence 
of the DOS over the whole bandwidth, as numerically 
evaluated from the detailed band dispersions, Eq. ([9]). 
In particular. Fig. U] shows p{E) for increasing strain, at 
fixed strain angle 9 — Q (armchair) and 9 = n/Q (zig zag). 
In both cases, for sufficiently low values of the strain 
modulus, the DOS depends linearly on E, according to 
Eq. (fT5|) . and the DOS slope increases with increasing 



strain, in agreement with Eq. (|19p and Fig. [3l However, 
while the spectrum remains gapless at all strains in the 
armchair case, a nonzero gap is formed at a critical strain 
in the zig zag case 9 = 7r/6, corresponding to the break- 
ing of the cone approximation at low energy. Such a 
behavior is confirmed by Fig. O showing the dependence 
of the DOS over the whole bandwidth, now at fixed strain 
modulus and varying strain angle. 

At sufficiently high energies, beyond the linear regime, 
the DOS exhibits Van Hove singularities both in the va- 
lence and in the conduction bands. As anticipated, these 
correspond to the occurrence of an ETT in the constant 
energy contours of either band dispersion relation -EkA, 
Eq. (I9|). As shown by Fig.lU the DOS is characterized by 
a single logarithmic cusp in each band in the unstrained 
limit (e = 0), that is readily resolved into two logarithmic 
spikes, both in the ^ = (armchair) and in the 9 — Tr/d) 
(zig zag) cases, as soon as the strain modulus becomes 
nonzero (e > 0). The low-energy spike disappears as 
soon as a gap is formed, corresponding to the breaking of 
the cone behavior around the Dirac point. Fig. [5] shows 
that the situation is indeed richer, in that the applica- 
tion of sufficiently intense strain along generic {i.e. non 
symmetry-privileged) directions allows the development 
of three logarithmic singularities in the DOS for each 
band, corresponding to the three inequivalent ETTs de- 
scribed in Section ini Again, the lowest energy Van Hove 
singularity disappears into the gap edge when the energy 
spectrum ceases to be linear around the Dirac points. 
This takes place when the Dirac points ik^) tend to ei- 
ther edge midpoint Mi of the IBZ {i ^ 1,2,3). In this 
case, the energy gap can be found explicitly and, in 
the simple case of no overlap {si = 0) , can be written as 



Ai = 2^tl + tl+tl^2tit2-2tit:i + 2t2h,{2Qa) 
As = 2^tl + tl+tl- 2tit2 + 2tit3 - 2t2i3,(20b) 
A3 = 2^tl +tl+tl + 2tit2 - 2tit3 - 2t2h, (20c) 

for i = 1,2, 3, respectively. 

Further insight into the anisotropical character of the 
low-energy cone dispersion relations around the Dirac 
points, Eq. p5|) can be obtained by recasting them in 
polar coordinates {q, (f), where q — {qcoscj), gsini/)). One 
finds therefore Eqx — v\{(j))q, the anisotropic prefactor 
vx{(j)) depending on the Dirac point around which one is 
actually performing the expansion. Fig. [6] shows wa(</') 
for the conduction band (A = 2) centered around ho- 
One notices that applied strain increases the anisotropy 
of the (j) dependence, until a critical value is reached, at 
which the cone approximation breaks down. This corre- 
sponds to a nonlinear behavior of Eqx along a specific 
direction (j)Q, characterized by the vanishing of v\{(j)) and 
given explicitly by 



C0tg(/)o 



3 ti±t2' 



(21) 
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FIG. 4: (Color online) Energy dependence of the DOS over 
the whole bandwidth, for increasing strain modulus e = ~ 
0.25 and fixed strain direction 6 = (top panel) and 9 = n/6 
(bottom panel). In both cases, the DOS slope close to the 
Fermi energy increases as a function of strain. However, while 
the DOS remains gapless for ^ = 0, a nonzero gap opens 
around £ = at a critical strain for 6 = 7r/6. 



when jial = \ti =F in Eqs. ([13]), and to the opening of a 
finite gap around zero energy in the DOS. In that case, 
the Fermi velocity vanishes along a direction (jjQ given by 



cotg^o = 



(l + £ll)cotg(/)o-£i2 
(1 + £22) -~ £2lCOtg0o' 



(22) 



IV. OPTICAL CONDUCTIVITY 

The paramagnetic component of the density current 
vector in momentum space reads^^^ 



dp 

(2^ 



(2p + p')c|,Cp+p', 



(23) 



where Cp (cj,) are destruction (creation) operators in the 
plane wave representation. In the homogeneous Hmit 
(zero transferred momentum, p' = 0), one has^ 
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FIG. 5: (Color online) Energy dependence of the DOS over 
the full bandwidth, for fixed strain modulus e = 0.25 and 
varying strain direction. 



where H is the system's Hamiltonian. 

Within linear response theory, the conductivity a is re- 
lated to the current-current correlation function through 
a Kubo formula 



huNAcell 



nL (0,0,0.), (25) 



where n is the electron density, T is the temperature, 
uj is the frequency of the external electric field, Aceii 
is the area of a primitive cell, and Ylf;^{p,p' ,uj) is the 
{I, m) component of the Fourier transform of the retarded 
current-current correlation tensor. One is usually con- 
cerned with the dissipative part of the conductivity ten- 



> 




jV(0) = _[i7,r] = -er, 



ih 



(24) 



V((|))/V°( COS(l 



FIG. 6: (Color online) Polar plot of v\{4>) (with A = 2, i.e. 
for the conduction band) around kjj, v\{(f)), normalized with 
respect to its value in the absence of strain, v'^p . Strain is here 
applied at a generic fixed angle 6 = 7r/4. The anisotropy of 
the Fermi velocity increases with increasing strain, until the 
shape of v\{(j}) breaks at £ = 0.28. This corresponds to the 
existence of a direction (solid blue line), Eq. (|2ip . along which 
the dispersion relation Eqx displays a nonlinear character. 
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sor, i.e. its real part. One has therefore and {p, r) denotes the Fourier transform of the param- 

agnetic component of the current density vector, at the 

aim{tJ.,T;uj) ^ -T— J-j Imn}J„(0,0,w), (26) imaginary time r. Projecting J^{p,t) onto the tight- 

ru-oNAceii binding states, and neglecting the Drude peak, one finds 

where is the retarded version of 

nim(p,p'.T) = -(r,[J7(p,T)j,^(p',o)l), (27) 



kA 



kA 



(28) 



where ^kA = -EkA ^ M ^^^d np{£^) denotes the Fermi function at temperature T. In the direction of the external field, 
i.e. for I — one finds 



M'^) 2Tr smh{^h(3\uj\) 1 ^ 2 

= — ;^ Z^|J; (k)|;^xi^(/3,/i;k)(5(fiw- (i^kA -£^kA)), 



Co 



(29) 



kA 



r 



where /3 = (fceT)"^ is the inverse temperature, J;^(k) = 
e^j;^(k), do = 7re^/(2/i) is proportional to the quantum 
of conductivity, Tq"^ = 16t^/(3A/37r?i^), and 



i^(/3,/i;k) 
^k 



2e'5(^'<-^)nj.(a,i)nF(Ck,2), (30a) 
n ' (30b) 



with = when (?k = (no overlap). 

We have numerically evaluated the longitudinal optical 
conductivity aii(uj), Eq. as a function of frequency 
w > at fixed temperature /cbT — 0.025 eV, for several 
strain moduli e and directions 9, as well as field orien- 
tations, here parametrized by the angle (j) between the 
applied electric field and the lattice x direction. Figs. [7] 
and [HI show our results in the case of strain applied in 
the armchair direction {9 ~ 0) . Fig. [7] shows a contour 
plot of the longitudinal optical conductivity au/ao as a 
function of frequency lo (radial coordinate) and applied 
field angle (polar angle). In the relaxed limit (e = 0), 
o'ii/ctq is isotropic with respect to the applied field angle, 
and exhibits a maximum at a frequency that can be re- 
lated to the single Van Hove singularity in the DOS (cf. 
Fig. [4]). Such a maximum is immediately split into dis- 
tinct maxima, in general, as soon as the strain modulus 
£ becomes nonzero. This can be interpreted in terms of 
applied strain partly removing the degeneracy among the 
inequivalent underlying ETTs. Such an effect is however 
dependent on the field direction 0, as is shown already by 
the anisotropic pattern developed by aujoQ in Fig. [71 for 
e ^ 0. Indeed, Fig. [5] shows plots of an jai^ as a function 
of frequency for fixed strain modulus e = 0.1 and vary- 
ing field orientation 4> = — 7r/2. The relative weight 
of the three maxima depends on the relative orientation 
between strain and applied field. Here and below, we 



consider the case = 0. A nonzero value of the chem- 
ical potential would result in a vanishing conductivity 
below a cutoff at uj k, smeared by finite temperature 
effects^. 



An analogous behavior is recovered when strain is ap- 
plied along the zig zag direction 9 = ir/Q, as shown in 
Figs. [9] and [101 Again, applied strain breaks down the 
original isotropy of the optical conductivity with respect 
to the field orientation in the relaxed case, with two max- 
ima appearing as a function of frequency (Fig. [9]). The 
optical weight of the different maxima depend in general 
by the relative orientation between strain and applied 
field. While the presence of the two peaks can be traced 
back to the existence of inequivalent ETTs, whose degen- 
eracy is here removed by applied strain, the last panel in 
Fig. [HI shows that at a sufficiently large strain modulus 
(here, e = 0.275), a gap opens in the low-energy sector 
of the spectrum, which is signalled here by a vanishing 
optical conductivity (dark spot at the origin in last panel 
of Fig. [SI . 



Finally, Figs. [TT] and [T2l show the longitudinal optical 
conductivity in the case of increasing strain applied along 
a generic direction, viz. 9 = 7r/4. Like in the previous 
cases, applied strain removes the isotropy of an j with 
respect to the field orientation (/). However, the degen- 
eracy among the three inequivalent ETTs is here lifted 
completely, and three peaks in general appear in the lon- 
gitudinal optical conductivity as a function of frequency, 
as shown also by Fig. [121 The redistribution of optical 
weight among the three peaks is now more complicated, 
as it in general depends on both the strain direction 9 
and the field orientation 6. 
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FIG. 7: (Color online) Polar plots of longitudinal optical con- 
ductivity aii/ao, Eq. (I29|l . as a function of frequency ijj > 
(polar axis) and electric field orientation </!> (azymuthal direc- 
tion). Here, we set fi = and ksT — 0.025 eV. Strain is 
applied along the 6 = (armchair) direction, and the strain 
modulus increases from left to right, and from top to bottom 
(e = 0,0.075,0.175,0.275). 



FIG. 9: (Color online) Same as Fig. [T] but for strain applied 
the 6* = 7r/6 direction. 




10 12 14 16 18 
ffl[eV] 

FIG. 8: (Color online) Longitudinal optical conductivity 
aii/ao, Eq. (|29|l . as a function of frequency uj > 0, for fixed 
strain modulus e = 0.1 and strain direction 8 = (armchair). 
Different lines refer to various orientations of the electric field 
(0 = 0,^/4,^/2). 



V. CONCLUSIONS 

We have discussed the strain dependence of the band 
structure, and derived the strain and field dependence 
of the optical conductivity of graphene under uniaxial 
strain. Within a tight-binding model, including strain- 
dependent nearest neighbour hoppings and orbital over- 
laps, we have interpreted the evolution of the band dis- 




e=0.1,(|)=0 

E=0.1,l])=3l/4 
E=0.1,ll)=3l/2 



2 4 6 



10 12 14 16 U 
a)[eV] 



FIG. 10: (Color online) Same as Fig. (8] but for strain applied 
the = 77/6 direction. 



persion relations with strain modulus and direction in 
terms of the proximity to several electronic topological 
transitions (ETT). These correspond to the change of 
topology of the Fermi line as a function of strain. In 
the case of graphene, one may distinguish among three 
distinct ETTs. We also recover the evolution of the loca- 
tion of the Dirac points, which move away from the two 
inequivalent symmetric points K and K' as a function 
of strain. For sufficiently small strain modulus, however, 
one may still linearly expand the band dispersion rela- 
tions around the new Dirac points, thereby recovering a 
cone approximation, but now with elliptical sections at 
constant energy, as a result of the strain-induced defor- 
mation. This may be interpreted in terms of robustness 
of the peculiar quantum state characterizing the electron 
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FIG. 11: (Color online) Same as Fig.[71 but for strain applied 
the S = 7r/4 direction. 




FIG. 12: (Color online) Same as Fig.[8l but for strain applied 
the 8 — TV /A direction. 



liquid in graphene, and can be described as an instance 
of 'quantum protectorate'. For increasing strain, two in- 
equivalent Dirac points may merge into one, which usu- 
ally occurs at either midpoint Mg = 1, 2, 3) of the first 
Brillouin zone boundary, depending on the strain direc- 
tion. This corresponds to the breaking down of hnearity 
of the band dispersions along a given direction through 
the Dirac points, the emergence of low-energy quasipar- 
ticles with an anisotropic massive low-energy spectrum, 
and the opening of a gap in the energy spectrum. Be- 
sides, we confirm that such an event depends not only 
on the strain modulus, but characteristically also on the 
strain direction. In particular, no gap opens when strain 
is applied along the armchair direction. 

We derived the energy dependence of the density of 



states, and recovered a linear dependence at low en- 
ergy within the cone approximation, albeit modified by a 
renormalized strain-dependent slope. In particular, such 
a slope has been shown to increase with increasing strain 
modulus, regardless of the strain direction, thus suggest- 
ing that applied strain may obtain a steeper DOS in the 
linear regime, thereby helping in increasing the carrier 
concentration of strained samples, e.g. by an applied 
gate voltage. We have also calculated the DOS beyond 
the Dirac cone approximation. As is generic for two- 
dimensional systems, the proximity to ETTs gives rise 
to (possibly degenerate) Van Hove singularities in the 
density of states, appearing as logarithmic peaks in the 
DOS.' 

Finally, we generalized our previous results for the op- 
tical conductivity^^, to the case of strained graphene. We 
studied the frequency dependence of the longitudinal op- 
tical conductivity as a function of strain modulus and 
direction, as well as of field orientation. Our main re- 
sults are that (a) logarithmic peaks appear in the optical 
conductivity at sufficiently high frequency, and can be 
related to the ETTs in the electronic spectrum under 
strain, and depending on the strain direction; (b) the 
relative weight of the peaks in general depends on the 
strain direction and field orientation, and contributes to 
the generally anisotropic pattern of the optical conduc- 
tivity as a function of field orientation; (c) the opening 
of a band gap, where allowed, is signalled by a vanish- 
ing optical conductivity. Thus, an experimental study of 
the optical conductivity in the visible range of frequen- 
cies as a function of strain modulus and direction, as 
well as of field orientation, should enable one to identify 
the occurrence of the three distinct ETTs predicted for 
graphene^. 

In our study, we have assumed that the chemical poten- 
tial docs not itself depend on strain. In the doped case, 
this is clearly an approximation, as the carrier concentra- 
tion is expected to remain constant, while the band struc- 
ture is modified by strain. It will therefore of interest, 
for future studies, to investigate the strain-dependence 
of the chemical potential required to maintain a constant 
carrier concentration. This will enable one to evaluate 
the dependence of the Hall resistivity on uniaxial strain, 
which is a quantity of experimental interes t^^i*^^ . 



APPENDIX A: SECTION OF A DIRAC CONE 
FOR STRAINED GRAPHENE 



Eq. P?|) yields the canonical form of the ellipse ob- 
tained as a section with constant energy E of the cone 
approximating the band dispersions around either Dirac 
point kc Eq. (jlSp. The center {qxQ,<lyo) with respect to 
k/) of the ellipse evolves linearly with energy E according 
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-0.025 -0.02 -0.015 -0.01 -0.005 0.005 0.01 0.015 0.02 0.025 

q,'°l a/E 



FIG. 13: (Color online) Showing the scaled position of the 
center (q^o /E, qyo / E) of an elliptical section of the Dirac cone 
around ko at constant energy E, Eq. (|A1|I . Each line refers 
to a given strain angle 6 with respect to the lattice x axis, 
and varying strain modulus e = — 0.2. E.g., the directions 
9 = and 9 = tt/G correspond to the vertical bottom and 
horizontal left line, respectively. 



to 



qxo = -A^ {d^ COS ri - dy sin ri)E 



(Ala) 



%o = - B"^ {d^ sin rj + dy COST]) E. (Alb) 



The ellipse semiaxes A, B are given by 
1 1 



1 

52 



(7- v/"'+/5') 



(7+ v/"'+/3')- 



(A2a) 
(A2b) 



In the above equations, we have made use of the foUowing 
definitions: 



COST] 



smr] 




2cosry\/a2 + 



7 



3\/3, 



a 



3^2 

— {tj+tl + tl), 



(A3a) 

(A3b) 

(A3c) 
(A3d) 
(A3e) 



where d is given in Eq. pop . One finds a,/? — > 0, while 
7 %t^a? /2 in the hmit of no strain, e 0. Fig. [T^ 
shows the dependence on the strain modulus of the scaled 
coordinates of the ellipse center for different strain orien- 
tations. 
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